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Abstract. We present a non-perturbative canonical analysis of the D = 3 quadratic- 
curvature, yet ghost-free, model to exemplify a novel, "constraint bifurcation", effect. 
\ Consequences include a jump in excitation count: a linearized level gauge variable 

i-Sh • is promoted to a dynamical one in the full theory. We illustrate these results with 

their concrete perturbative counterparts. They are of course mutually consistent, as 
CO " are perturbative findings in related models. A geometrical interpretation in terms of 

■ propagating torsion reveals the model's relation to an (improved) version of Einstcin- 

\ Weyl gravity at the linearized level. Finally, we list some necessary conditions for 

CO " triggering the bifurcation phenomenon in general interacting gauge systems. 
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1. Introduction 

Canonical analysis a la Dirac is a straightforward, if sometimes labyrinthine, approach 
to counting a system's physical degrees of freedom (DoF) in the presence of gauge 
symmetries and non-linear interactions. However, this approach can uncover unexpected 
subtleties, as already exemplified by some toy models in [1] . In this paper, we show that 
more physically motivated theories can also contain similar subtleties, with qualitatively 
important consequences. Our focus will be on a specific self-interacting spin-2 gravity 
model, but similar effects may well arise in other interacting theories with higher-spin 
gauge symmetries, at least if some (listed) necessary conditions are met. 

The theory we shall study in detail is the truncation of D — 3 "NMG" [2] to its 
pure quadratic curvature, yet ghost-free, part [3], 

I[g] = ^J d 3 x^ G^G^ - l -G 2 =^J d 3 xy^ G, U S^, S, u : = R^—g^R 

(1) 

here G^ v is the Einstein tensor, G is its trace and S^ v is the D = 3 Schouten tensor. 
We have set an overall dimensional constant to unity and used mostly plus signature. 
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This model provides a remarkable example of symmetry-breaking through the 
clash between its two local, conformal and coordinate, invariances. Their co- 
existence at linearized level underlies this fourth derivative metric system's "miraculous" 
transmutation into single ghost-free vector excitation, or equivalently to a propagating 
torsion with non-propagating metric. Nonlinearly, however, conformal- is necessarily 
sacrificed to coordinate-invariance. 

This paper is organized as follows: In section [2] we perform the Hamiltonian analysis. 
In section [3] we exhibit the bifurcation mechanism and count the number of physical 
degrees of freedom. In section HI we first transmute the free field into a geometric 
model with propagating torsion, relating it to an (improved) version of Einstein- Weyl 
gravity, as well as to its other, Maxwell vector, avatar. We then exhibit the nonlinear 
obstructions and their effects, in particular, introduction of propagator-less variables. 
In section [5] we comment on possible applications to other interacting gauge theories. 



2. Hamiltonian analysis of Schouten gravity 

In keeping with our Hamiltonian approach, we use a first-order formulation of the action, 



I[e,u,f, A] = J d 3 x 



(2) 



Roman/Greek indices are local/world. In form notation, the Cartan, dreibein e* (with 
determinant e) and (dualized) spin-connection u\ variables define the torsion and 
(dualized) curvature T l = de l +£ l jk oo^e h , R l = duj 1 + e % jkU! : 'uj k . The Lagrange-multiplier 
A* ensures the on-shell torsion constraint T l = 0, while the auxiliary f l is essentially the 
same as in [3] and in the Appendix A of the second reference [I]. The first-order 
action (T2]) is just an Ostrogradsky auxiliary variable form of the original fourth-order 
action ([1]), so the two are classically equivalent and share the same excitation content. 

We now proceed to analyze ([2]) canonically, following the earlier methods of 
[HE], to which we refer for a more extensive discussion. The Lagrangian's variables 
(e* M , u\, X\, f\) have conjugate momenta (tt^, U.^, Pi u ', P^), leading to the primary 
constraints (~ means weakly equal, i.e., equal on the constraint surface): 

0,° := tt, w , K := - e 0a/3 A l/3 w , (3a) 

:= IT w , := IL a - e 0a/ % w , (3b) 

Pi"«0, Pf w . (3c) 

The constraints (<fri a , §i a ,pi a , Pi a ) are second class. It will be useful to define the linear 
combination 

k° ■■= + + w . (4) 

Eliminating the momenta (7Tj Q , IL ,^ , Pj Q ) leads to a partly reduced phase space, in 
which the nontrivial Dirac brackets are given by 

{e\, AV = rf'eoafi , K Q , f j p} = r/ ij e 0Q/3 . (5) 
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The remaining Dirac brackets are the same as the corresponding Poisson brackets. The 
canonical Hamiltonian T-L c can be written as a sum over secondary constraints 

H e := e\Ui + uj\fC t + f\K t + A 4 7- (6) 

up to boundary terms. Here, 

Hi ■= -e 0a ^V a X^ + e L°V K - \ U(f k0 - fe k0 )^j « , 



(7a) 



Ki := -e 0a/3 (V a f ip + e ijk e i a \ k f l ) « , (7b) 



and 



Hi ■= -\p*Ru# + \ e (fi° ~ fei°) w , (7c) 
7I:=-^e°^T fa ^0, (7d) 

Vk := J - / 2 ) • (8) 

The consistency conditions of the secondary constraints lead to the ternary constraints 
6(iu ~ ~ V'/iK that establish the symmetry of the auxiliary fields. Similarly, the 
consistency conditions of the ternary constraints yield quaternary constraints 

X = A^0, ^ = /+^«0. (9) 

At this stage the constraint procedure fortunately ends: no further constraints are 
generated through consistency conditions. 



3. Counting degrees of freedom: Bifurcation 

Having found all constraints in the previous section, we establish their first /second class 
properties in order to count the number of physical DoF. We summarize the main result 
in tables [T] and HJ The fact that we have two different tables is a consequence of the 
advertised subtlety we called "bifurcation" and now explain. 

The consistency condition of the quaternary constraint ip leads to the following 
Dirac bracket 

{<p, H T } = - A CT - fgT) z»u w , (10) 
where 1-Lt is the total Hamiltonian 

U T := H c + u'o^i + v\^ G + w\ Pl ° + z\P? (11) 

and (u, v, w, z) are Lagrange-multipliers of primary constraints. All components 
and z^y are determined at this stage (and hence the corresponding constraints are second 
class), except for the Lagrange multiplier zqq that multiplies the primary constraint 
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P 00 . The consistency conditions fllOp allows us to determine the remaining Lagrange 
multiplier z 00 , unless the condition 

/ 00 = fg 00 (12) 

holds. The resulting bifurcation consists in the following: if condition (fT2l) does/not 
hold, the quaternary constraint (p and the primary constraint P 00 are both first/second 
class. 





First class 


Second class 


Primary 


A , 


„.0 p.O 


Secondary 






Ternary 




#0/3, #a/3, "00^9, 


Quaternary 







Table 1. Classification of constraints in the partly reduced phase space absent the 
constraint (TT^j) . 





First class 


Second class 


Primary 


A , P 00 


Pi , remaining Pj° 


Secondary 






Ternary 




#0/3, #c>!/9, "^0/3, ^a/3 


Quaternary 




X 



Table 2. Classification of constraints in the partly reduced phase space with the 
constraint ((T^|) present. 



According to tabled! i.e., when ( 1T2|) is absent, we have a 48-dimensional phase space 
with 12 first class and 20 second class constraints. Consequently, the theory in general 
exhibits two local physical DoFs, namely the massive bulk gravitons. However, if the 
condition ( I12jl holds, then table |2] states that there are 14 first class constraints and 18 
second class constraints. In that case the number of local physical DoFs is reduced to 
one. 

Translating condition (1T2]) into metric form, by using f^ u = 2S^ IU and / = |P, one 
finds that it holds automatically if the metric is a solution of the Einstein equations 

Rt » _ gl»> R = Q ^ f^ =g ^f_ ( 13 ) 

This observation concurs with the fact that the linearized theory around Einstein 
solutions has one additional gauge symmetry due to partial masslessness [6]. Again, 
this symmetry enhancement results from the fact that the constraint P 00 becomes first 
class, leading to the additional gauge symmetry. However, the latter is an artifact of 
linearization, being broken in the full nonlinear theory, as we have shown above. 
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4. Conformal versus coordinate invariance 

In this section we provide an alternative derivation of the loss of gauge-invariance beyond 
linearization in a straightforward perturbative approach. This sheds additional light 
on the bifurcation mechanism and exhibits it as a clash between diffeomorphism and 
conformal invariance. Moreover, the perturbative analysis can be useful particularly for 
possible higher-rank/higher-spin generalizations, where a full canonical analysis is often 
less accessible than a perturbative one. 

We work with the second order, "Ostrogradski" action [equivalent to (pQ)], using 
auxiliary, symmetric tensor density, variables f^ u : 

I[9,f] = \J d'*{GMr-l[&-&f)^/V=g} (14) 

where we have omitted the contracting metrics. [Completing squares and integrating 
out / recovers ([I]).] We will study ( Tl4|) . initially at linear, then full non-linear, metric 
levels. The linearization of ([T]) is manifestly (linear) diffeo-invariant, while use of the 
Bianchi identity easily confirms its conformal invariance, under 

Shftv = -2r)fj, v a, SG^ = {d^dy - g^ u O)a , (15) 

with our convention = d\T x ul j, — d u T x \ /1 + . . . . Both invariances also hold in (114"]) . of 
course, with / transforming as a (linear) diffeo tensor and conformally like the Schouten 
tensor: Sf^ = d^d u a. 

Now we count DoF: Varying h^ u yields G^ y {f) = 0, where G is the usual linear 
Einstein operator. In D = 3, there are no Einstein excitations (Riemann and Ricci 
being equivalent), so f^ u is a pure gauge "metric". Varying /, we learn that h^ u obeys 
the Einstein equation with linear source, so the general solution is: 

/am; = dpAv + duA/t , GftuQh) = f,,,, - 7]^f ■ (16) 

Inserting - legally - (fTBT) into (fT4"]l . we note first that its G(h) f = 2G^ u d^A„ term 
vanishes by the Bianchi identities upon part integration, leaving the quadratic /-terms: 
these precisely combine into the promised reduced one-DoF Maxwell action (second 
reference in [2]): 

I[h, /] -> ~\ J d 3 x F, u := (d,A u - djjAfj) ; (17) 

it is invariant under SA^ = \d^a. That a is indeed our conformal transformation 
parameter 

8 fit* = d^dva (18) 

is then verified by (I15II16I) . The above set of field equations is consistent with both 
underlying invariances; for example, upon taking the divergence of Einstein equation in 
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(ITB|) . its left side vanishes by the Bianchi identity, while the divergence/ conservation of 
its (symmetric) right side matter source also does, being proportional to the latters's 
field, i.e., Maxwell's, equations. In this connection, note that there is no "spin- loss" 
paradox in the above tensor-to- vector transmutation because all massless fields in D = 3 
are necessarily spinless [HIE], at least in flat space; it would be useful to learn if this 
persists in (A)dS. 

The above, conformal-to-gauge transmutation should not be confused with a 
separate, surprising [9], D = 3 conformal invariance enjoyed by Maxwell, by virtue of its 
further transmutability into a scalar. We provide a concise derivation, emphasizing the 
nonlocal (as usual with such transmutations) price involved: The first order Maxwell 
action is 

I[F, A] = - 1 -Jd 3 x [F^{d,A y - d v A») - \Fl v \ , (19) 

where F^ u and are independent variables. Varying A^ gives d v F^ u = 0, whose 
general solution is F^ u = e fl ' /a d a S, the scalar S having dimension of A^; F's dual vector 
*F I1 is just dyS, a relation whose highly non-local inverse we also note, 

S(*F) = cr 1 ^ * F^ . (20) 

Replacing *F by dS in ( !T9|) immediately yields the free scalar action, I[F, A] — > 
I[S] = —\ J d 3 x(dS) 2 . To cancel the scalar stress-tensor's trace, add to its minimal 
T^S) the usual identically conserved (hence not affecting the Poincare generators) 
improvement [10] term, ST^ ~ —\{d^d v — rj^DjS 2 . Its on-shell Tr(dS) 2 cancels that 
of T^y. The above procedure is nonlocal in F^ u , inherited from S(*F) in (1201 . But one 
cannot improve directly at the vector level either. Any such attempt is doomed from 
the start: it is easily verified that there is no identically conserved 5T ~ dd(AA) - as 
required by dimension - that is even on-shell gauge- invariant, so one cannot cancel the 
Maxwell stress-tensor's trace, ~ \F 2 . Nor is any other local choice available, since it 
must, just by dimensions, depend on AA rather than on the field strengths FF; the 
scalar itself is the only gauge- invariant, dimensionally correct, but nonlocal, option. [We 
have not attempted to find a direct descent from ( |T4l) to the scalar action.] 

In canonical terms, the associated DoF-reducing Maxwell gauge constraint is the 
standard (A V- E), leaving & D = 2 + 1 photon with just one transverse DoF, whose 
sign is fixed by that chosen for (CQ) or ([HI). More explicitly, the Maxwell action's kinetic 
term is (— E T -A T — E L -A L ), in terms of the spatial transverse-longitudinal orthogonal 
vector decomposition Vi = e^djV + diV, which commutes with time-derivatives. The 
longitudinal excitation is removed by the Gauss constraint (A V ■ E) that enforces 
E L = 0. Instead, the cubic correction ( |23~1) contains terms quadratic in Aq, thereby 
replacing this constraint by an irrelevant (because it integrates out) perfect square plus 
(cubic) terms that depend on (E\Ai), so the longitudinal DoF are reinstated. (We 
have checked that cubic terms ~ hoo(A ) 2 are absent though.) Not having studied the 
Hamiltonian in detail, we cannot assert that it is no longer bounded below, but that 
seems likely for any cubic: one would have to include quartic corrections for a meaningful 
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conclusion, though that is not very relevant any more. This paradoxical second order 
form of a fourth order action is explained by the original metric propagator indeed 
being [3] ~ (V 2 D) -1 , thus agreeing with / having dimension of curvature, / ~ ddh, 
hence A ~ dh. [Actually, our method also applies, with appropriate dimensional 
numerical differences, to the D = 4 Weyl action, which can also be written (modulo 
coefficient differences in Schouten) in the form (I14p . The D = 3 argument is unchanged, 
in particular that h and F each obey the Einstein equation, so each has 2, rather than 
D = 3's 0, DoF. The A^, coordinate vector gauge, part of / still represents a (now 2 
DoF) photon [3] (see also [IT]). Again, there is no spin-transmutation paradox (despite 
the presence of spin D — 4), because the propagator of a fourth derivative order theory 
can (and usually does) include lower-spin poles.] 

Our linearized system has a second, more general geometrical interpretation: it 
represents a propagating vector torsion and contortion, but a non-propagating (because 
D — 3) metric. Indeed, this is already suggested by the Einstein equation in (TIE]) . 
We recall that one may extend Riemannian geometries by extending its metric affinity 
r%i/(/i) to include an antisymmetric torsion, T"^, and a contortion, whose symmetric 
part we denote by K a ^ v ) [we define (anti-)symmetrization with a factor 1/2]. For our 
system, the choices are 

TV = b S^A V] , KV) = c S^A v) + d ^ v A a . (21) 

The parameter choice d + 2c = kills the Ricci tensor's antisymmetric part. The further 
choices b = 3 and c = 1 make the vanishing of the Einstein tensor, but now of the full 
connection, G(T tot ) = 0, T tot a pu := r%(/i) + T a ^ u (A) + K a {pu) (A), coincide with the 
source-ful Einstein equation of (1T61) . The divergence of the right hand side of that 
equation is proportional to the Maxwell operator, whose vanishing is then assured by 
consistency with the metric Bianchi identity. Note that this total connection r to t is 
metric compatible in the Einstein-Weyl sense. 

{V li {h)-2A ll )g af ) = (22) 

In this interpretation the Maxwell field is recognized as the Weyl potential, up to a 
factor and the presence of torsion, see e.g. [12]. We stress that the inclusion of torsion 
leads to an improvement of the usual torsionless Einstein-Weyl connection, since our 
Ricci tensor is symmetric, while the Einstein-Weyl Ricci-tensor in general is not. More 
explicitly, for Einstein-Weyl we have b = and c = —2d; note that metric-compatibility 
in the Einstein-Weyl sense is guaranteed if b + c + 2c? = 0. Interestingly, the tracefree 
part of the right Eq. ( TT6|) coincides precisely with the Einstein-Weyl equations (up 
to terms quadratic in the Weyl potential), see p2J[13] and Refs. therein. This provides 
yet-another-way to see why conformal invariance is broken cubically: the Einstein-Weyl 
equations are Weyl-invariant, but they coincide with our equations only at the linearized 
level, not non-linearly. 

We now show explicitly how enforcing diffeo- destroys conformal- invariance beyond 
linear order, as it must since the full action (CQ) involves the factor (^—gg^g 013 ) ~ 
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0{g~ 1 / 2 ), rather than O(g ). We need only consider the first, cubic, deviation, where 
the effect will be manifested as loss of Maxwell gauge invariance. Returning to ()14p and 
inserting the linearized values ( 1T6|) . we find that the cubic action reduces, schematically, 
to the "bare A^-form" 

J 3 [/V, d " x l 2G * ih)d»A v + h^Q^(A)] , (23) 

where (and also its integral) is an (irreducibly) gauge variant quantity bi-linear 
in the gauge potential A^ This manifest loss of Maxwell gauge invariance means 
of course, loss of conformal invariance. A concrete, "no-calculation" realization of 
invariance loss is now easy. Consider first the vacuum state in the gauge g^ v = n^, 
A^ = 0: all linearized and cubic terms manifestly vanish. Now gauge-vary this "null" 
configuration, with local parameter a = x^x^, where x M are the Minkowski coordinates. 
Then A^ = x M , f^ u = f = 6, h^ v oc x 2 i]^ u : This means A^ is a conformal Killing 

vector of the flat background. The (gauge invariant) quadratic action still vanishes of 
course. Instead, J 3 oc J t] f j 1 uQ flu - This, manifestly constant trace, rj^Q^, is easily found 
to be non-zero; hence the cubic action is a non-invariant, proportional (in this gauge) 
to the vector "mass" term J 3 oc = m 2 A fJ- A fl ^ 0, already at vacuum. 

We conclude that the simple one-DoF content of the linearized level is indeed lost 
here, with the implied consequent presence of additional, ghost, modes that always 
plagues generic quadratic actions, albeit without propagators of their own. 

5. Conclusions 

We have studied a novel, higher-order, clash between two local invariances that 
characterize dynamical quadratic curvature Schouten Schouten-ghost-free D = 3 gravity, 
at linear order. After performing a canonical analysis and exhibiting a bifurcation 
in the constraint analysis that leads to different counts of DoFs, we considered the 
theory perturbatively, providing a transparent derivation of the free theory's one-DoF, 
second order, character: our geometrical, torsionful, representation complementing an 
earlier vectorial transmutation. We then traced the unavoidable breaking, already at 
cubic level, of the model's conformal symmetry by its nonlinear diffeo- invariant, but 
conformal-gauge dependent, completion. The culprits were the dynamical (rather than 
Minkowski) metrics that contract indices beyond lowest order. 

The above symmetry-breaking raises one, also novel, field-theoretic puzzle posed by 
this otherwise consistent model: Its propagator depends on fewer variables than do its 
vertices. How does one calculate (at least perturbatively) either classically or at quantum 
loop level? The new variables, having no lines of their own, can only lie on open, but 
presumably not on closed loop lines, yet they are not external fields either. We have also 
checked the consistency of our perturbative analysis with a non-perturbative canonical 
analysis along the lines of jl], with perfect agreement. A (vaguely) similar situation 

1 We obtain = A a df,(f^r) a P - rTf aP ) - W(.f a pf afi + I 2 ) + \fF v - f» a F a v . 
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occurs in topologically massive gravity, whose metric's propagator is not uniformly of 
either second or third derivative order, since its third derivative Cotton contribution is 
independent of the metric's conformal factor, one that is present in its Einstein, second 
order term. That problem is not one of principle, however - all components have proper 
propagators, just ones of different momentum order - and can only affect topologically 
massive gravity's UV behavior p3]. Rather, the nearest analog is perhaps massive D = 4 
Einstein gravity with a "wrong" explicit mass term a la [15]; however that symmetry- 
breaking pathology is inserted by hand, rather than, as in our model, from enforcing a 
greater invariance! Clearly, some intriguing unsolved directions remain. 

The general non-perturbative conditions of the canonical analysis we encountered 
here are a non-perturbative extension and here agree with the perturbative results. 
Moreover, a similar effect was observed in generalized massive gravity, see second 
reference in [I]. This suggests that the bifurcation effect we described here is not 
just confined to our specific model ([1]), but can be present in more general interacting 
higher-rank/ higher- spin gauge theories as well as, perhaps, in vector models, see [16] 
and Refs. therein. 

We conclude with a list of ingredients necessary, but not always sufficient, for the 
emergence of bifurcation. 

(i) Higher (at least spin 2) gauge symmetries 

(ii) Non-linear (at least cubic) interactions 

(iii) The linearized theory must have additional (linearized) gauge symmetries such as 
conformal /Weyl symmetry 

(iv) The linearized gauge symmetry must be broken in the interacting theory 

It would be interesting to construct explicit examples of bifurcation in D > 3 theories 
with spins higher than 2, and/or towers thereof. 
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Note added: As this work was completed, and simultaneous with our submission of the 
perturbative analysis (essentially section H]), the closely related posting [17] appeared. 
It discusses the Chern-Simons and other extensions of the present model. Where they 
overlap, our conclusions agree; indeed the models of [17] provide further evidence for 
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our conjecture that the bifurcation we encountered here is a more generic feature of 
interacting higher-rank/higher-spin theories. 
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